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ABSTRACT 

For extrasolar planets with orbital periods, P>10 days, radial velocity surveys find non-circular 
orbital eccentricities are common, (e) ^ 0.3. Future surveys for extrasolar planets using the transit 
technique will also have sensitivity to detect these longer period planets. Orbital eccentricity affects 
the detection of extrasolar planets using the transit technique in two opposing ways: an enhancement 
in the probability for the planet to transit near pericenter and a reduction in the detectability of 
the transit due to a shorter transit duration. For an eccentricity distribution matching the currently 
known extrasolar planets with P>10 day, the probability for the planet to transit is ~ 1.25 times 
higher than the equivalent circular orbit and the average transit duration is ^ 0.88 times shorter than 
the equivalent circular orbit. These two opposing effects nearly cancel for an idealized field transit 
survey with independent photometric measurements that are dominated by Poisson noise. The net 
effect is a modest ~ 4% increase in the transiting planet yield compared to assuming all planets have 
circular orbits. When intrinsic variability of the star or correlated photometric measurements are the 
dominant source of noise, the transit detectability is independent of the transit duration. In this case 
the transit yield is ~25% higher than that predicted under the assumption of circular orbits. Since the 
Kepler search for Earth-sized planets in the habitable zone of a Solar-type star is limited by intrinsic 
variability, the Kepler mission is expected to have a ~25% higher planet yield than that predicted 
for circular orbits if the Earth-sized planets have an orbital eccentricity distribution similar to the 
currently known Jupiter-mass planets. 

Subject headings: eclipses — planetary systems — techniques: photometric 



1. INTRODUCTION 

The known extrasolar planets possess a b road distri- 
bution of orbital eccentricity f|Butler et al.l [2006) (see 
Figure [1]). The short period. Hot Jupiter (P<10 day) 
planets predominately have circular orbits. However, at 
longer orbital periods circular orbits become a minority 
and the median eccentricity for extrasolar planets e~0.3. 
At the extre me ecce ntricity end, there are t hree planets, 
HP 8 0606b (jNaef e t al. 200 l|), HP 2078 2b (jJones et all 
120061 ). and HP 4113b (,Tamuz et al.ll2007[ ). having e>0.9. 
HP 80606b comes closer to its stellar host (a=0.033 AU) 
than many o f the circular orbit Hot Jupiter planets. 
iFord fc Rasiol (j2007[ ) recently reviewed the various mech- 
anisms invoked to explain the distribution of eccentrici- 
ties for the known extrasolar planets. Interactions with 
a stellar companion, planetary companion, passing star, 
gaseous disk, planetesimal disk, and stellar jets have all 
been proposed to modify the orbital eccentricity of ex- 
trasolar planets. 

Limited discussions in the literature have been given to 
the impact orbital eccentricity has on a trans it survey for 
extrasolar planets. iTinglev fc SacketH (|2005f l discuss the 
impact of eccentricity on their rj parameter (77 is the ra- 
tio of the observed transit duration to an estimate of the 
transit duration). As expected, they find transits occur 
near pericenter (apocenter) are shorter (longer) in dura- 
tion than the circular orbit case, and they show that the 
transit duration of an eccentric is typically shorter than 
a circular orbit of the same period. However, their dis- 
cussion was focused o n the impact of orbita l eccentricity 
on their 77 parameter. iMoutou et al.l (|2006f ) also discuss 
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how the transit duration is affected by orbital eccentric- 
ity, but they do not quantify the irn p act th is will have 
on transit surveys. Recently, iBarnesI (|2007t l derives the 
probability for a planet on an eccentric orbit to transit 
and conclude that the photometric precision of current 
surveys and future surveys, such as Kepler, is insuffi- 
cient to determi ne the o rbital eccentricity solely from 
the light curve. iBarnesI (|2007t ) concludes that without 
knowledge of the eccentricity from radial velocity data 
or independent measurement of the stellar host radius, 
the habitability of planets detected with Kepler will re- 
main unknown. 

Neglecting the impact eccentricity has on transit de- 
tections is justified for the current sample of transit- 
ing planets given the predominance of circular orbits 
for the Hot Jupiters and the strong bias of transit sur- 
veys against finding lon g period planets on circular or- 
bits (|Gaudi et al.ll2005l ). The announced transit for the 
planet orbitin g HP 17156^ with a 21 day period and ec- 
centric orbit ( Barbieri et al.l 120071 ) is a precursor for the 
kinds of planets detectable in transit surveys. As transit 
surveys continue, longer period transiting planets may 
be discovered. More importantly, the recently launched 
COROT mission will surpass current surveys for sen- 
sitivity to longer period planets (P~ several months; 
iBorde et al.ll2003l ). Also, the Kepler mission, scheduled 
for launch in 2009, has a goal to find Earth-sized ob - 
jects at 1 AU from their host star (jBorucki et al.l[2b04( l. 
The main purpose of this paper is to show that given the 
distribution of eccentricities for the currently known ex- 

^ F irst detected by the radial velocity technique l|Fischer et al.l 
[20071) 
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trasolar planets, eccentricity should not be ignored in as- 
sessing the detectability of transiting giant planets when 
the transit survey is sensitive to planets with P> 10 day. 

In addition to longer period planets, the COROT 
and Kepler missions also will detect transiting plan- 
ets with small, Earth-sized radii. The eccentricity 
distribution for planets less massive than the cur- 
rently known J upiter-mass planets is beginning to 
be explored iRibas fc Miralda-Escudel (|2007) and 
iFord fc: Rasiol (|2007 ^ provide tentative evidence for the 
tendency of lower mass planets to have lower eccentric- 
ities in the current sample of radial velocity planets. 
There is theoretical agreement that a proto-planetary gas 
disk strongly damps the eccentr icity of non-gap open- 
ing e mbedded low-mass planets (|Goldreich fc Tremaind 
119801 ICresswell et al.. ,2007) . However, the theoretical 
models show that once the gas disk dissipates, dynamical 
interactions amongst the planets results in a random- 
walk diffusion that leads to an increasing eccentricity 
that takes on a Rayleigh-distribution simil ar to what 
is ob s erved (dotted line in Figure [TJ Juric fc Tremaind 
I2007t IZhou et all [20071) . Opposing the mcreases m ec- 
centricity from planet-planet scattering, late stage inter- 
actions with planetesimals can preferentially damp the 
eccentricities of lower mass planets enabling theoretical 
models to achieve the low ecc entricities of the terrestrial 
planets of the Solar System (|Ravmond et al.l 2006) and 
possibly explain the cur rent trend of lower e ccentricities 
for lower mass planets (jFord fc Rasiol [2007h . Given the 
number of potential physical processes that can affect 
orbital eccentricity, it is premature to assume that the 
typical Earth-like planet has an eccentricity near zero 
like the Solar System. 

In this study, §[2] reviews the observed eccentricity dis- 
tribution of the known radial velocity extrasolar planets. 
The broad distribution of orbital eccentricity has two 
main effects on the sensitivity of a transit survey. First, 
the planet-host separation varies along an eccentric orbit, 
enhancing the probability to transit when the planet is 
relatively closer to the stellar host. §[3] quantifies the net 
affect on the transit probability for a population of plan- 
ets with non-circular orbits. Second, the planet velocity 
varies along and eccentric orbit, resulting in a reduction 
or lengthening of the transit duration. § |4] quantifies the 
distribution of transit durations resulting from a popula- 
tion of planets on eccentric orbits. § [5] describes how the 
transit duration affects transit detection for transit sur- 
veys in the limit of various noise sources. § [6] concludes 
by quantifying the net result of the two aforementioned 
effects, the enhanced probability to transit and the re- 
duced detectability, on the yield from transit surveys. 

2. ECCENTRICITY DISTRIBUTION 

The solid histograms in Figure [1] show the normal- 
ized distribution of eccentricities fo r the known extra- 
solar planets with F > 10 day (Butl er et aLl |200tf). The 
top panel shows the eccentricity distribution from plan- 
ets reported before November 2006 and the bottom panel 
shows the eccentricity distribution including more recent 
discoveries (September 2007). The more recent radial 
velocity discoveries over the last year have increased the 
sample of low eccentricity systems. Results from this 
study are given for both epochs of the observed eccen- 
tricity distribution in order to characterize how the un- 
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Fig. 1. — Normalized histogram of the eccentricity distribution 
for known extrasolar planets with P>10 day. The eccentricity dis- 
tribution is given for the known extrasolar planets at two epochs: 
November 2006 {Top Panel) and September 2007 {Bottom Panel). 
For this study a two-piece model parameterizes the eccentricity 
distribution {solid line). The model has 2 parameters; Ccrit and 
Gmax- The model distribution is uniform at low eccentricity up to 
Ecrit; and toward higher eccentricity, the model linearly decreases 
to zero- value at fimax ■ Also shown in the top panel is the Rayleigh 
distribution that describes the eccentricity distribution from the 
planet-planet scattering calculations of Quric & Tremaine 200'3) 
{dotted line). The dotted histogram in the bottom panel shows the 
eccentricity distribution from the latest epoch (September 2007) 
that excludes planets with P< 20 day. 

certainty in the underlying orbital eccentricity distribu- 
tion affects the results. No attempt was made to select 
planets with well determined eccentricities, avoidance 
of multiple planet or multiple star systems, correction 
for observationa l biases against high eccentricity planets 
(|Cumming||200l . or evolution of orbital elements due to 
tidal circularization. Planets with P < 10 day are not 
included in order to examine the potential issues that 
future transit surveys sensitive to longer period planets 
will encounter. In this period regime (P>10 day) high 
eccentricity planets are quite common ((e) = 0.3 and 
ae = 0.2). 

A two-piece model parameterizes the eccentricity dis- 
tribution from the earlier epoch (top panel of Figure [1]). 
The first piece at low eccentricity is flat up to Ccrit- The 
second piece matches the first piece at ecrit, linearly de- 
creases toward higher eccentricities, and is zero at emax- 
The normalized two-piece model is given by the equation. 



I Ecrit +emax " _ ^ii., 

I Cmax < e < 1 

(1) 

A minimization yields the best model parameters 
Gcrit — 0.25 and Cmax = 0.92. The minimization 
applied to the more recent epoch eccentricity distribu- 
tion (bottom panel of Figure [1]), yields ecrit = 0.0 and 
Smax = 0.91. The relative increase in the number of low 
eccentricity planets discovered in the last year results in 
b est param eters that are ef fectiv ely a single-p iece model. 

lJuric fc T remaind (|2007| ) and IZhou et al.l ([2007) pre- 
dict that the eccentricity distribution of dynami- 
cally active planetary systems approaches a Rayleigh- 
distribution due to planet-planet scattering, and they 
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show tlie Rayleigli-distribution is similar to the observed 
planet eccentricity distribution. The dotted curve in Fig- 
ure [1] shows the Rayleigh-distribution with ae = 0.3 
that best describes the out come of the plane t -plane t 
scattering calculations from lJuric fc Tremain^ (|2007l ). 
The Rayleigh-distribution under represents the low ec- 
centricity systems in the most recent eccentricity data. 
lJuric fc Tremainj (|2007D based their study on the sam- 
ple of planets known at the time (April 2006). Subse- 
quently, relatively more planets with low eccentricities 
have been announced. In addition, lJuric fc Tremaini 
()2007h included systems with P>20 day whereas this 
sample has P>10 day. The dotted histogram in the lower 
panel of Figure [1] shows the eccentricity distribution at 
the most recent epoch (September 2007), but with sys- 
tems with P<20 day removed; apparently the difference 
is slight and will hereafter be ignored. Although this 
study concentrates on the two-piece model to describe 
the eccentricity distribution, results are also given for 
the Rayleigh-distribution with ae = 0.3. To provide re- 
sults with the Rayleigh-distribution, the distribution is 
set to zero for e > 0.95. 

3. TRANSIT PROBABILITY 

iBarned ()2007| ) derives the impact of orbital eccentric- 
ity on the pro bability for a planet t o transit its stellar 
host (see also ISeagroves et al.l l2003| ) . The probability 
to transit depends on the planet-star separation during 
transit. For an eccentric orbit, a transit occurring during 
pericenter enhances the transit probability and a transit 
occurring during apocenter decreases the transit prob- 
ability. When averaged over observing angles, the net 
result is an enhancement of the probability to transit 
over the circular orbit case. 



Probxe = 



Probxo 
(l-e2)' 



(2) 



where Pro bxf) is the tra nsit probability of the circular 
orbit case ()Barnesl 120071 ) . A planet with e=0.6 is ^--^1.5 
times more likely to transit than a planet on a circular 
orbit. 

With the observed eccentricity distribution from § [51 
it is possible to calculate the average enhancement in the 
transit probability over the circular orbit case. 



(Probi 



ProbTeP(e)de. 



(3) 



The resulting integral using Equation [T] for the eccentric- 
ity distribution is. 



(ProbTc) = 



PiobT 



[2(e„ 



Ccrit )arctanh(eci:it ) 



-I- In 



- In 



/ (l-ecrit)(l 
V(l-Hecrit)(l 



(4) 



Using the numerical values for Ccrit and e,nax from 
§ [21 (ProbTe)=l-26 and 1.23 for the earlier and cur- 
rent epochs of observed eccentricities, respectively. The 
Rayleigh distribution results in (ProbTe)=l-31. Figure[2l 
illustrates (Probxe) for other choices of the eccentricity 
distribution model parameters. Figure [2l is a function of 
Cmax and the curves are for selected values of ecrit as la- 
beled. Thus, it is expected that the yield from a transit 
survey that is sensitive to P>10 day planets will be ^^25% 
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Fig. 2. — Average transit probability as a function of the observed 
eccentricity distribution model parameters. The transit probability 
is scaled to the equivalent transit probability for a circular orbit. 
The abscissa indicates Emax , and the curves are for selected values 
of eciit as labeled. For the model parameters shown in Figure [T] 
(Prob-pe) ~ 25% higher than assuming all planets have circular 
orbits {dotted lines). 

larger than expectations that assume circular orbits for 
all planets. This section discusses only the probability 
to transit, and it does not address whether the planet 
has a transit signal that is detectable. Thus, the re- 
sults from this section assumes the enhanced probability 
for the planet to transit does not affect its detectability, 
which is the topic of the next section. 

4. TRANSIT DURATION 

4.1. Edge-on Transit 

In addition to enhancing the probability for a transit 
to occur, orbital eccentricity results in the transit dura- 
tion varying according to the planet's longitude of peri- 
center during transit. The transit duration is shortest 
(longest) if the transit occurs when the planet is at peri- 
center (apocenter) when transiting its stellar host. A 
short transit duration may reduce the detectability of the 
transit event. Assuming constant velocity during tran- 
sit, the transit duration scaled to the edge-on (i — 90°) 
circular orbit case has extrema of 



ts. 



'(ire) 

(lie) 



iTe- 



72, 



(5) 



where the sign on top and bottom corresponds to the 
ericenter and apocenter transit durations, respectively 
Barnes! [20071) . When e—0.6, the transit duration at 
apocenter is twice the circular orbit case. Given the 
range of eccentricities observed for extrasolar planets 
with P>10 day, we expect significant variations in the 
tr ansit duration compared to that of a circular orbit. 

iTinglev fc SackettI ()2005f ) (their Equation 7) provides 
a simplified form of the transit duration, including lon- 
gitude of pericenter, nj, and orbital inclination, In 
deriving the transit duration, (jTinglev fc Sacke"ttil2005| ) 
assume constant orbital velocity and planet-star separa- 
tion during transit and the planet crosses the stellar disk 
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resulting in 
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Fig. 3. — Left: Edge-on transit duration, t, scaled to the edge-on 
circular orbit case as a function of the longitude of pericenter, ro, 
for several eccentricities as labeled. Right: For a uniform distribu- 
tion of ro, solid curves show the distribution for transit duration 
scaled to the edge-on circular orbit case assuming all orbits are 
edge-on (i.e. i = 90°). 



along a straight path (a ^ R^). An exact calculation 
(of a numerical nature given the need to solve Kepler's 
Equation) of the transit duration finds Equation 7 of 
iTinglev fc Sacketti (|2005,) is accurate to better than 5% 
for e=0.9 and orbital separation a > 0.1 AU (P~10 day). 
The equation is accurate to better than 1% for e=0.9 and 
a > 1.0 AU. 

After scaling to the transit duration of the edge-on (i = 
90°), circular orbit with the same perio d and separation, 
Equation 7 of iTinglev k SackettI (|2005D becor 



3mes, 



(1 + e cos(ti7)) 



l-p2 



(l-e2) 



(1 + ecos(ti7)) 



cos(i)2 



(6) 

where p = a/ {R* + Rp), R* is the radius of the star and 
Rp is the radius of the planet. The definition for zu is 
with respect to the line of sight. Thus, w = means 
the pericenter is aligned with the observers line of sight, 
and zu — 180° means the apocenter is aligned with the 
observers line of sight. This varies from the definition 
of the argument of pericenter, to, which is defined with 
respect to the line of nodes on the plane of the sky {-uj = 
w- 90°). 

The i — 90° case illustrates the first order impact of or- 
bital eccentricity on transit duration. In this case, Equa- 
tion [HI simplifies to 



'90" 



(1 -I- ecos(n7)) 



(7) 



The left panel of Figure [3] shows the transit duration 
with respect to the edge-on circular orbit case for a va- 
riety of orbital eccentricities as a function of zu. For a 
uniform distribution in ro, the shallow slope in transit 
duration at pericenter and apocenter results in a dis- 
tribution of transit durations peaked at these extrema 
(given by Equation [S]) with low probability for a transit 
duration in between these extrema. 

Assuming a uniform distribution of zu, p{zu)dzu = 
dzu/iT, the transformatio n law of probabilities (Equation 
7.2.4 in lPress et al.|[l99^ yields the distribution of tran- 
sit durations. 



p{T)dT 



dzu 



p{zu)dT, 



(8) 



p{T)dT = 



VT^dT 



(9) 



The right panel of Figure [3] shows the probability density 
for transit duration with i — 90° for various orbital ec- 
centricities. As expected the probability density is heav- 
ily weighted toward the extrema. The singularities in the 
probability density at the extrema are integrablc, making 
the probability density normalizable. 

4.2. Affect of Varying Inclination Angle 

The previous section shows the expected transit du- 
ration for fixed orbital eccentricity when the inclination 
angle, i = 90°. Varying the inclination impacts the tran- 
sit duration in two important ways. First, as the inclina- 
tion decreases, the path of the planet across the stellar 
disc shortens, resulting in a reduction of the transit dura- 
tion. Second, for an eccentric orbit, the planet is closer 
to the star at pericenter and farther away at apocen- 
ter. The smaller separation at pericenter increases the 
range of inclination angles that result in a transit, and 
conversely, a planet at apocenter will have a reduction 
in the range of inclination angles for a transit to occur. 
Thus, long transit duration events when the planet is 
at apocenter become increasingly rare as the eccentric- 
ity increases. This section quantifies these effects on the 
transit duration distribution. 

Beginning with a simple joint distribution that is uni- 
form in cos(i) and uniform in vj enables the more com- 
plicated distribution for transit duration to be derived. 
For this section, e remains fixed. The beginning joint 
distribution is, 

P{zu,{) < cos(i) < cos{imin)\e-)d{cos{i))dzu — Ad{cos{i))dzu, 

(10) 

where cos(zmm) is cosine of the minimum inclination an- 
gle necessary for a transit to occur, and A is the normal- 
ization constant. Setting r = in Equation[5]and solving 
for cos(i) yields cos{imin) = (1 + e cos(n7))/p(l — e^). In- 
tegrating Equation [10] over the range of variables yields 
the normalization constant, A = p{l — e^)/7r. 

The transformation law of probabilities for multiple 
dimensions (Equation 7.2.4 in lPress et al.|[l99^ provides 
the joint probability density in terms of new variables r 
and zj when starting with the probability density for 
cos(i) and zuf. 



P{t, zu\e)dTdzu 



cScos(z) cS cos(z) 



dr 



P(cos(i), zu/)dTdw, 



(11) 

In Equation [Til dzu/dr = and dzu/ /dzu = 1. Thus, 
only dcos{i)/dT remains, and Equation 1 1 1 1 simplifies to 



P(t, zu\e)dTdzu 



9cos(i) 



p{l~ e'^)/7TdTdzu. (12) 



The requisite derivative needed in Equation [12] is ob- 
tained by solving Equation [6] for cos(z) yielding 



cos(i) = 



(1 + ecos(ci7)) 
p(l-e2)3/2 



V(l-e2)-r2(l + ecos(n7))2, 

(13) 
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as tlie solid line. The function, 




Fig. 4. — Left: Transit duration distribution scaled to the edge- 
on circular orbit case for a uniform distribution of cos(j) and fixed 
eccentricity. The cases e=0.2, 0.4, 0.6, and 0.8 are shown. Right: 
Average transit duration scaled to the average transit duration of 
a circular orbit as a function of orbital eccentricity {solid line). 

and taking the partial derivative with respect to r. Per- 
forming this operation yields 

P(0 < r < Tq, TUmhi < < Tr\e)dTd'cu 
r(l -t- ecos(n7))^ 



ttVI - e2^(l - e2) - r2(l -f- e cos(n7))' 



zdrdvo, 



(14) 



where the lower limit to tu is necessary to avoid an imag- 
inary result and is given by 



MIN 1.0, 



vr 



(15) 



For T < Tperi, ^min = 0, but aS T > Tperi , ^min > 

since only transits occurring closer to apocenter result in 
a transit duration long enough. 

Finally, the conditional probability density for transit 
duration alone is obtained by integrating over 



P{T\e)dT 



P{t, ■uj\e)dTdu:!. 



(16) 



A solution to the above integral is possible in terms of 
a summation of the incomplet e elliptic integrals of th e 
first, second, and third kind (|Bvrd fc Friedman! Il954f ). 
In practice, I choose to solve the integral numerically, 
which is readily solved using the Romberg open ended 
algorithm which takes into account the si ngularity at 
the lo wer limit of integration at tu — tUmin ( Press et al.l 
[1991) . Tests of convergence show Equation [T5] has a sin- 
gularity at r = Tpej.j, but it is integrable elsewhere over 
the range < r < Tq. 

The left panel in Figure |4] shows the probability den- 
sity for transit duration scaled to the edge-on circular 
orbit case at several values of eccentricity. By compari- 
son to Figure [31 the main impact of orbital inclination is 
to strongly enhance the probability of observing a short 
duration event at pericenter relative to a long duration 
event at apocenter. The probability density also allows 
arbitrarily short events due to the potential for grazing 
events. 

The probability density for transit duration can be 
summarized by finding the average transit duration at 
fixed eccentricity. 



[ TP{T\e)dT. 
Jq 



(17) 



The right panel in Figure [Jj shows (te) scaled to the aver- 
age transit duration of the circular orbit case, (tq) = 7r/4, 



1^ 



(18) 



fits the relation to better than 10~^ (similar to the nu- 
merical integration precision), which very strongly sug- 
gests this is t he analytical solut i on to the integral in 
Equation [TTj iTinglev fc Sacketd (|2005f ) find the same 
result in terms of the average value for their r] parameter 
at fixed eccentricity (see their Equation 18). 

4.3. Transit Duration Distribution For Observed 
Eccentricity Distribution 

The previous section derives the transit duration dis- 
tribution at fixed eccentricity. This section derives the 
transit duration distribution assuming planets follow the 
observed eccentricity distribution. As in the previous 
section, deriving the transit duration distribution begins 
with the simple distribution that is uniform in cos(i), uni- 
form in w, and the distribution of e follows the observed 
distribution as given in Sj2j The transformation law of 
probabilities enables transforming the simple distribu- 
tion in cos(i), wf, and e/ into the distribution expressed 
in r, VJ, and e. 

The initial distribution is given by, 

P(0 < cos(i) < cos(i™„), < n7/ < TT, < e/ < e^nax) 
dcos{i)dwidei — AP{ei)d{cos{i))dwidei , 



(19) 



where the normalization constant A = /t^, where 



7 = (e 



In 



2 

max 
(1- 



[2{e 



(1-eLt) 



+ e„ 



In 



ecrit)arctanh(ecrit) 
(l-ecrit)(l + e 



(l + e„.it)(l 



(20) 



The Jacobian transformation matrix simplifies as before, 
such that the transit duration distribution scaled to the 
edge-on, circular orbit transit duration is given by. 



P(0<T<T„n7 

d cos(i) 



dr 



mill — ^ } ^min ^ ^ ^ ^max 

P(cos(i), n7/, ei)dTdwde, 



)dTdvDde 



(21) 



where the lower limit to the eccentricity, Cmim becomes 
necessary for r > 1, when too small of an eccentricity 
cannot produce a transit duration as long as r. Solving Ta 
for e yields e^in = MAX[0.0, (r^ - 1)/{t^ + 1)]. Overall, 
the joint distribution is given by 

P(r, vj, e)dTdvDde — 
7P(e) t(1 + ecos(cc7))^ 



(1 - e2)3/2^i e2 - r2(l + e cos(tu))2 



zdrdvude. 
(22) 



Integrating over vj and e provides the final probabil- 
ity density for r for the assumed distribution of orbital 
eccentricities. Given the additional complication of in- 
tegration over two variables, an analytical solution was 
not forthcoming. As in § 14.21 the singularities in the 
integrand arc integrable, and the Romberg open ended 
algorithm which takes into account the singularity at the 



6 



Christopher J. Burke 



1.5 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



A 

V 0.80 
0.75 
0.70 



0.0 










0.6 \\\ 




0.7<VV\- 




0.8 \V 




0.85\\ 




0.90 



^max 



Fig. 5. — Left: Transit duration distribution scaled to the edge- 
on circular orbit case for a uniform distribution of cos(j) and the 
observed distribution of e as shown in the top panel (solid line) 
and bottom panel (long dashed line) of Figure [T] Given the bias 
against detecting high eccentricity planets in radial velocity sur- 
veys, a uniform distribution of orbital eccentricities up to high 
eccentricity (e^rit = 0.9 and emax = 0.95) will be highly skewed 
towards t ~ 0.3 (short dashed line). A Rayleigh distribution of 
orbital eccentricities results in relatively fewer transit durations 
T ~ 1 due to the relatively fewer objects on circular orbits (dot- 
ted line). Right: Average transit duration scaled to the average 
transit duration of a circular orbit as a function of the observed 
eccentricity distribution model parameters. The abscissa indicates 
Gmax, and the curves are for selected values of ecrit as labeled. For 
the model parameters shown in Figure [T] (r) ~ 0.88 times shorter 
than assuming all planets have circular orbits (dotted line). 

lower hmits of integration (jPress et al .119921 ) provides the 
sohition. 

The sohd and long-dashed lines in Figure [5] shows the 
distribution of transit duration scaled to the edge-on cir- 
cular orbit case for a population of extrasolar planets 
that follows the observed orbital eccentricity distribu- 
tions shown in top and bottom panels of Figure[Tl respec- 
tively. The probability density has a broad flat top with 
r = 0.4 just as likely to occur as r = 1.0. For compar- 
ison, the transit duration distribution that would result 
from a nearly uniform distribution of orbital eccentricity 
up to high eccentricities (ccrit = 0.9 and e,„ax = 0.95) 
is shown as the short-dashed line in Figure [HI Given the 
bias against detecting e > 0.6 extrasolar planets in radial 
velocity studies (Gumming 2004), a large population of 
e ^ 0.9 planets cannot be ruled out. The transit dura- 
tion distribution resulting from the Rayleigh distribution 
of orbital eccentricity (dotted line in Figure [5]) has rela- 
tively fewer planets with t ^ 1 due to its relatively fewer 
circular orbits. 

In the future, if a statistically large sample of transiting 
planets with orbital P > 10 days is available with accu- 
rate stellar parameters, histograms of observed transit 
duration scaled to the edge-on circular orbit case may 
help characterize the underlying eccentricity distribu- 
tion. After accounting for the selection effects, a large 
number of r '--^ 0.3 detections relative to r ~ 0.8, as il- 
lustrated in the right panel of Figure [5l would indicate 
e = 0.9 planets are as common as circular orbits. Work 
toward understanding the sensitivity of Kepler for con- 
straining the underlying eccentricity distribution is un- 
derway (E. Ford, private communication). 

The right panel of Figure [5] summarizes the transit 
duration distribution by showing (rg) in terms of the pa- 
rameters characterizing the eccentricity distribution, Ccrit 
and Cniax- Each line corresponds to a fixed value of ecrit 
as labeled, and the abscissa indicates emax of the eccen- 
tricity distribution. For the eccentricity distribution in 



the top panel of Figure [T] {t^) / (to) — 0.88 and the bot- 
tom panel {te) / {tq) = 0.89. The Rayleigh distribution 
results in {Te)/{To) — 0.86 

5. DISCUSSION: APPLICATION TO TRANSIT SURVEYS 

The results from § [4] quantify the impact orbital ec- 
centricity has on the transit duration. Overall, orbital 
eccentricity results in shorter transit durations than the 
circular orbit case, and the short transit duration reduces 
the transit detectability. This section quantifies the re- 
duction in transit detectability for various noise models 
of transit surveys. 

In a transit survey with independent photometric mea- 
surements, the transit signal to noise ratio is, 

SNR- y]^, (23) 

a 

where AF is the transit depth (the transit is modeled as a 
box-car shape), a is the error on a photometric measure- 
ment, and iVobs is the number of measurements during 
one or more transit (s). A shorter transit duration re- 
duces iVobs oc T, thus, SNR cx ^/T. The observed eccen- 
tricity distribution (Figure [T]) results in f = {Te) / {tq) — 
0.88, and on average the reduced transit duration results 
in an effectively smaller SNReS = Vt — 0.94 per transit 
than assuming all planets are on circular orbits. 

The above impact on the transit signal SNR due to 
a shorter transit duration is for an individual star in a 
survey. However, the reduced SNR will have a larger 
impact on the overall transit detec tability in a n idea l 
field transit surve y. As described in iGaudi et al.l ()2005[ ) 
and iGaudil ()2007| ). a specified SNR criteria for transit 
detection, SNRmin, in a field transit survey corresponds 
to a maximum distance, ^max oc SNRj^^jj^, out to which a 
planet is detectable. This proportionality assumes white 
noise and the dominant source of photometric err or is 
Poisson noise. In the studies of ' Gaudi et al.l (|.2005) and 
iGaudil (|2007f l. ^max is a function of the planet radius and 
stellar host spectral type (i.e. fmax is a smaller distance 
for a smaller radius planet or larger radius star) . For this 
study only the dependence of £max on transit duration is 
of interest. 

Overall, the yield from a transit survey is propor- 
tional to the number of objects in the survey that meets 
SNRmin, which is TVobj oc i'^-^^^ for stars distributed uni- 
formly in the survey volume as appropriate for nearby 
stars. The effective SNRmin^eff = SNRminV^ larger than 
assuming all planets are on circular orbits. Thus, the 
number of objects in an idealized transit survey where a 

transit is detectable is A^obj.e = '?^''^A^obj,o = 0.827Vobj,o 
times smaller than the case where the detectability of 
a transit is based on assuming all planets have circular 
orbits, iVobj.o- Despite the reduced detectability of tran- 
sits, this is offset by the higher probability for the planet 
to transit in the case of significant eccentricity § [3l The 
overall yield of the idealized transit survey is discussed in 
§ [6] taking both the reduced detectability and enhanced 
probability to transit into account. 

In practice, transit surv eys typically are affected by 
correlated measurements (jPont et al.l |2006[ ). In this 
regime, the correlation time scale is similar to the tran- 
sit duration and repeated measurements do not add 
independent information. When correlated measure- 
ments dominate the photometric error, the SNR = 
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(AF/crcor)\/^tr, wfiere iVtr is tlie number of transits de- 
tected, and the correlated measurement error, tTcor, no 
longer depen ds on the stellar luminosity, but is constant 
(|Gaudi|[2007f ). When dominated by correlated measure- 
ments, the SNR is independent of r and the only re- 
quirement is to have a short enough sampling cadence 
to detect the shortest transit duration expected. Thus, 
a shorter r due to orbital eccentricity has no impact on 
the transit detectability in a survey which is dominated 
by correlated measurements (i.e. iVobj.c = -^obj,o)- How- 
ever, the Poisson limited, white noise transit survey will 
have a higher planet yield than a survey d ominated by 
correlated measurements (jPont et al.ll2006l ). but in the 
latter case, orbital eccentricity does not impose any ad- 
ditional reduction in the transit detectability. 

The Kepler space-based transit su rvey whose goal is to 
detect several Earth-sized planets (jBorucki et al.l [200^ 
contends with stell ar intrinsic va riability as the domi- 
nant source of noise ()Jenkinsll2002l ) . Using integrated flux 
measu rements of the Sun to model intrinsic variability of 
stars, iJenkin^ ( |200^ characterizes the detectability of 
Earth-sized transiting planets with Kepler. The Sun has 
low noise on 7 hr time scales typical of transiting Earth- 
sized planets in 1 yr orbital periods. However, the solar 
noise increases by ~ 4 orders of magnitude by 10 day 
time scales (see Figure 2 of Jenkins 2002). The rapidly 
increasing intrinsic noise toward longer time scales can- 
cels any benefit of increased transit signal SNR for transit 
durations >4 hr for the Kepler mission (see Figure 9 of 
lJenkinsll2002l ). Thus, as long as the transit duration re- 
mains above 4 hr, orbital eccentricity will not reduce the 
detectability of Earth-sized planets with Kepler. 

6. CONCLUSION 

Orbital eccentricity results in an enhanced probabil- 
ity for a planet to transit and potentially a reduction 
in the transit detectability. The overall yield from a 
transit survey is given by A'dot Probx x A^obj- The 
results from this study can be used to scale the over- 
all yield from a transit survey based on assuming all 
planets are on circular orbits for an assumed distribu- 
tion of orbital eccentricity. The enhanced probability for 
a planet to transit, Probx with a distribution of orbital 
eccentricity scaled to the circular orbit case is given by 
Equation |4l The reduced number of transiting planets 
detectable A'obj scaled to the circular orbit case for an 
ideal transit survey where the photometric noise is white 
and dominated by Poisson error is given in § [H Multiply- 
ing these two factors provides the overall yield of an ideal 
transit survey scaled to assuming all planets are on circu- 
lar orbits. Figure [6] show that these two opposing effects 
nearly cancel over the parameters that characterize the 
eccentricity distribution. Thus, for an idealized transit 
survey with the currently observed orbital eccentricity 
distribution, the overall yield will be 4% greater than 
assuming all planets are on circular orbits. The result 
for an idealized transit survey with a Rayleigh distribu- 
tion of orbital eccentricities gives a similar enhancement 
(4%). 

In ground-based transit survey s, correlated mea sure- 
ments limit transit detectability (jPont et al.ll2006l ) and 
intrinsic variability of the star will limit the dete ctability 
for E arth-sized planets for the Kepler mission (jjenkinsl 
1200^ . In both cases the transit detectability is inde- 
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Fig. 6. — The overall yield from an idealized transit survey 
as function of the observe eccentricity distribution model parame- 
ters scaled to the assumption that all planets have circular orbits. 
The transit survey is assumed to have independent photometric 
measurements that are dominated by Poisson noise. The abscissa 
indicates emax, and the curves are for selected values of e^rit as 
labeled. For the model parameters shown in Figure [l] N^^i ~ 4% 
higher than assuming all planets have circular orbits {dotted line). 



pendent of the transit duration, in which case A^obj is 
independent of the orbital eccentricity. For these cases, 
the transit survey will have higher returns by a factor 
of (Probxe) (Equation [4|) than estimated by assuming 
all planets are on circular orbits. However, the reduced 
planet yield in a transit survey due to intrinsic stellar 
variability or correlated measurements must be properly 
accounted for. If every dwarf star has an Earth-sized 
planet orbiting in the habitable zone, then assuming cir- 
cular orbits, the Kepler mission expects to detect 100 
Earth-sized planets in the habitable zone dBorucki et al.l 
l2004f ). The work presented here indicates Kepler will 
have (Probxe) ~ 25% higher yield if Earth-sized planets 
in the habitable zone have a planet eccentricity distribu- 
tion similar to the currently known sample of giant plan- 
ets from radial velocity surveys. However, if Earth-sized 
planets with e~0.9 are as common as e^-^^O then, the yield 
of Earth-sized planets could be 80% higher from the Ke- 
pler mission assuming the high-e, short transit duration 
planets are still detectable. 

The dependence of the transit survey yield on the un- 
certain underlying orbital eccentricity distribution im- 
plies an uncertainty in measuring the frequency of ter- 
restrial planets in t he habitable zone ( a major goal of 
the Kepler mission iBorucki et al.l l2004f ) . An analysis 
of the transit yield from a transit survey that assumes 
all planets are on circular orbits will overestimate the 
frequency of habitable planets if high eccentricities are 
common and not taken into account. In practice a va- 
riety of noise regimes affect a transit survey and accu- 
rate yields necessitate an accurate understanding of the 
photometric noise, stell ar sample, and underlying eccen- 
tricity distribution (,Burke et al.ll2006l: iGould et al.1l2006t 
iFressin et aLll2007[ ). 
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